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ABSTRACT 

This  paper  presents  the  theory  and  examples  of  performance  for  a  new  algorithm  that 
tracks  targets  using  a  Multiple  Model  Unified  Joint  Probabilistic  Data  Association 
(MM-UJPDA)  filter.  The  models  in  the  MM-UJPDA  can  be  set  to  the  ambiguity  velocities 
encountered  when  initiating  tracks  on  a  sensor  that  has  ambiguous  velocities  in  its 
measurements.  Alternately,  the  models  can  be  set  for  tracking  manoeuvring  targets.  Thus 
each  parallel  filter  in  the  MM-UJPDAF  is  assigned  one  of  a  range  of  possible  target  model 
parameters.  The  term  'unified'  summarizes  a  number  of  key  features  in  the  algorithm. 
These  are:  multiple  non-uniform  clutter  regions,  a  model  for  a  visible  target  to  compute 
track  confidence  for  track  promotion,  and  measurement  selection  based  on  a  fixed  number 
of  nearest  measurements.  The  filter  formulation  used  a  new  approach  to  create  track 
clusters  for  determining  the  nearby  tracks  that  share  measurements.  The  filter's 
performance  is  demonstrated  with  track  initiation  using  the  multiple  model  and  multiple 
target  approach  while  for  established  tracking  only  the  multiple  target  approach  is  used. 
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Executive  Summary 

The  filter  described  in  this  report  extends  the  Multiple  Model  Unified  Probabilistic  Data 
Association  Filter  (MM-UPDAF)  that  is  currently  used  for  track  initiation  on  data  from 
the  Jindalee  Over-the-Horizon  Radar  at  Alice  Springs.  The  term  "Unified"  is  now  used 
to  summarize  the  following  features  that  are  included  in  this  tracker:  an  event  for  track 
initiation  and  termination;  the  selection  of  the  nearest  measurements  to  each  tracking 
filter's  predicted  position;  and,  an  a  priori  model  for  the  clutter  distribution  which  is 
represented  by  multiple  clutter  zones  that  have  both  uniform  and  nonuniform 
distributions. 

The  current  tracking  system  uses  the  MM-UPDAF  to  resolve  velocity  ambiguities  during 
track  initiation  and  when  a  track  is  promoted  to  the  established  state  the  UPDAF  is  used. 
The  extension  to  the  MM-UPDAF  that  is  reported  here  is  part  of  a  staged  development 
activity  to  further  improve  performance.  The  MM-UPDAF  and  the  UPDAF  do  not 
consider  the  possibility  that  measurements  can  be  from  another  target.  That  is,  they  treat 
tracks  as  isolated  entities.  Thus  in  dense  target  conditions,  instances  occur  where  a  track 
can  transfer  to  another  nearby  target.  Also  a  track  estimate  can  be  midway  between 
closely  spaced  multimode  returns  with  the  same  Doppler  and  range  but  different 
azimuth.  At  this  stage,  the  adaption  of  the  MM-UPDAF  to  handle  multiple  manoeuvre 
models  has  not  been  made.  This  was  not  done  because  targets  which  manoeuvre  are 
more  likely  to  do  so  in  the  proximity  of  other  nearby  targets.  Therefore  since  the 
MM-UPDAF  and  the  UPDAF  do  not  include  other  targets,  no  attempt  was  made  to 
enhance  manoeuvre  performance  over  the  current  capability. 

The  extension  to  the  MM-UPDAF  contained  in  this  paper  involves  a  new  target  model 
which  defines,  about  each  track,  a  cluster  of  nearby  tracks.  Thus  each  track  becomes  a 
reference  track  and  the  other  tracks  that  have  selected  the  same  measurements  are  linked 
to  form  a  track  cluster.  With  the  inclusion  of  multiple  target  models,  excessively  large 
track  clusters  are  formed  in  dense  target  conditions.  This  problem  lead  to  the  definition 
of  a  new  clustering  rule  to  constrain  the  growth  of  a  cluster  about  a  reference  track.  Even 
so,  cluster  sizes  increase  to  about  10  and  the  computation  of  the  number  of  measurement 
to  target  assignments  or  events  becomes  large.  To  evaluate  these  events,  a  joint  event  tree 
which  factorizes  the  events,  gives  an  efficient  implementation  which  allows  real-time 
tracking. 

The  common  name  given  for  a  filter  of  this  type  is  the  Multiple  Model  Joint  Probabilistic 
Data  Association  Filter  or  MM-JPDAF[8].  While  this  filter  has  already  been  formulated, 
the  approach  taken  here  significantly  differs  in  the  ways  described  above.  It  also 
incorporates  the  advanced  features  in  the  MM-UPDAF,  namely  the  measurement 
selection  approach,  the  inclusion  of  track  initiation  and  the  multiple  non-uniform  clutter 
zones.  So  in  this  paper,  this  new  multiple  target  filter  is  termed  the  MM-UJPDAF. 


This  filter  was  implemented  to  rim  on  the  Jindalee  replay  system  and  limited  testing  was 
undertaken  using  datasets  that  contained  closely  spaced  targets.  These  were  particularly 
challenging  to  the  current  tracker.  The  implementation  only  used  the  MM-UJPDAF  for 
track  initiation  while  for  promoted  tracks,  it  used  the  UJPDAF.  Thus  no  attempt  was 
made  to  try  multiple  manoeuvre  models  for  promoted  tracks.  At  this  stage  the  primary 
emphasis  has  been  to  improve  track  initiation  and  tracking  in  dense  target  conditions. 
Specific  examples  showed  how  this  new  filter  correctly  tracked  closely  spaced  targets. 
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1.  Introduction 

The  fundamental  problem  of  target  tracking  is  measurement  association.  When  a  sensor  does 
not  provide  any  identification  of  the  source  of  measurements,  the  measurement  source  must 
be  inferred  in  order  to  obtain  target  state  estimates.  When  the  environment  provides  a  high 
degree  of  clutter,  the  tracking  system  has  to  differentiate  between  target  measurements  and 
those  produced  by  interference  and  noise.  When  the  sensor  collects  measurements  from  several 
targets,  the  tracking  system  must  determine  which  measurements  were  caused  by  each  of  the 
targets. 

One  solution  to  the  problem  of  tracking  in  a  cluttered  environment  is  the  Probabilistic  Data 
Association  Filter  (PDAF)  [8],  The  PDAF  assigns  an  association  probability  for  each 
measurement  selected  by  a  track  and  uses  these  probabilities  to  weight  these  measurements  for 
track  update.  The  PDAF  allows  for  the  case  that  the  target  may  not  have  formed  any 
measurement.  This  is  due  to  either  a  failure  to  detect  the  target  signal  in  noise  or  the  target  not 
existing.  The  nonexistent  target  is  modelled  as  an  invisible  target.  Using  the  concept  of  target 
visibility  [1]  it  is  possible  to  formulate  a  PDAF  which  performs  automatic  track  initiation  and 
deletion. 

The  original  PDAF  formulation  has  some  limitations.  It  assumes  that  all  measurements  are 
caused  by  either  the  track  being  updated  or  by  clutter.  Under  dense  target  conditions,  this 
assumption  is  violated  and  the  effect  of  other  target  tracks  must  be  considered.  The  Joint  PDAF 
(JPDAF)  [2]  is  used  to  account  for  measurements  from  other  tracks. 

One  radar  that  forms  measurements  with  the  preceding  characteristics  is  Over-the-Horizon 
Radar  (OTHR).  Australia's  OTHR,  Jindalee  is  a  scanning  sensor  that  provides  measurements 
consisting  of  range,  azimuth  and  Doppler.  The  Doppler  measurement  can  be  ambiguous  due 
to  the  radar  operating  parameters.  There  are  also  objects  that  give  target-like  returns,  such  as 
transponders,  which  do  not  move  in  space  but  provide  an  observed  Doppler  shift.  To 
successfully  deal  with  these  ambiguities,  the  current  Jindalee  tracker  [1]  uses  a  PDAF  with 
target  visibility,  a  nonuniform  clutter  model  and  multiple  velocity  model  initiation. 

In  1998  a  JPDAF  was  implemented  within  the  Jindalee  Over  the  Horizon  Radar  (OTHR)  signal 
processing  replay  frame  work  [4].  This  filter  was  used  only  for  tracks  that  were  already 
established  and  the  current  system  [1]  was  used  for  track  initiation.  This  report  presents  an 
extension  of  this  JPDAF  to  cope  with  some  of  the  problems  of  track  initiation  noted  in  [1],  The 
JPDAF  extension  incorporates  the  multiple  model  initiation  scheme  of  [1]  as  well  as  the 
nonuniform  clutter  model.  The  aim  of  this  extension  is  to  improve  track  initiation  performance 
by  accounting  for  measurements  caused  by  nearby  tracks. 

This  area  of  initiation  on  interfering  targets  occurs  with  closely  resolved  multi-path  propagation 
for  OTHR.  In  this  case,  similar  propagation  paths  are  observed  as  parallel  tracks  and  may  be 
poorly  tracked  by  a  PDAF.  It  is  important  to  correctly  separate  these  multi-path  tracks  because 
conversion  from  the  radar  measurement  space  into  geographic  co-ordinates  requires  a 
knowledge  of  the  ionospheric  path  for  that  track. 
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2.  Outline  of  Target,  Sensor  and  Update  Model 


The  proposed  multiple-model  JPDAF  with  target  visibility  (MM-UJPDAF)  is  intended  to  be  an 
evolution  of  the  current  1RSU  multiple  model  PDAF  with  target  visibility  (MM-UPDAF). 
Therefore,  the  model  used  for  the  target  and  sensor  are  identical  as  that  used  in  the  existing 
MM-UPDAF  tracker  [1],  This  model  is  repeated  in  this  document  for  the  purposes  of 
completeness  and  clarity. 


The  model  used  to  represent  the 
environment,  target  and  sensor  is  shown  in 
Figure  1.  Here  the  sensor  collects  returns 
from  targets  and  the  environment  via  the 
propagation  medium.  These  returns  are 
converted  into  sensor  measurements  in  the 
receiving  and  measurement  processes. 
The  sensor  measurements  are  estimates  of 
all  target  like  signals  received  by  the  sensor 
for  each  time  step  k.  To  limit  the  number  of 
measurements,  a  threshold  is  used  to  pass 
likely  target  measurements.  Here  the  list  of 
measurements  is  assumed  to  be  dominated 
by  measurements  from  clutter  and  noise. 
The  source  of  a  measurement,  ie  target  or 
clutter,  is  also  assumed  to  be  unknown. 
The  measurement  process  also  determines 
the  extent  of  multiple  clutter  regions  and 
updates  the  clutter  model  parameters  that 
are  used  later  in  the  track  update  process. 


Figure  1.  Sensor  and  environment  model 


The  selection  process  in  Figure  1,  uses  the  track  parameters  in  the  track  data  base  to  select  those 
measurements  that  are  likely  to  be  from  the  target.  All  targets  that  give  measurements  in  the 
received  sensor  data  at  time  k  are  assumed  to  be  represented  by  tracks  in  the  track  data  base. 
The  contents  of  this  data  base  are  updated  by  the  measurements  and  the  clutter  model 
parameters  from  time  k  to  form  new  estimates  that  are  applied  to  the  data  at  the  next  time  k+ 1. 

2.1.  A  Priori  Data 

Before  outlining  the  processing  steps  in  Figure  1,  we  define  the  a  priori  models  for  the  target, 
the  clutter  and  track  data  base. 

2.1.1.  Target  Dynamic's  Model 

The  target  is  assumed  to  be  represented  by  a  fixed  number  of  M  discrete  target  models  that 
describe  the  target  dynamics.  The  mlh  target  model  for  target  track  t  is  given  by. 


Xn,,t  <*  +  1>  =  Fm,fM  X,nJ  (k)  +  U  (k),  m  =  1, . .  ,M 


(1) 
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Where  xm  t(j)  is  the  state  of  target  track  t  with  model  m  at  time  /,  F m  t(k)  is  the  state  transition 
matrix  for  the  m A  target  model  of  track  t,  and  uml(k)  is  the  mth  target  model's  process  noise 
vector,  Gaussian  distributed  with  zero  mean  and  covariance  Q(  m(k).  The  process  noise  accounts 
for  target  manoeuvres  and  model  errors.  Here  F mt(k)  and  Q m  t(k)  are  assumed  to  be  known  for 
each  model  of  each  track. 

The  target  model  switching  is  a  Markov  process  represented  by  the  a  priori  MxM  transition 
probability  matrix  with  transition  probabilities  plm„  where  /  is  the  model  index  at  time  k- 1  and 
m  is  the  model  index  at  time  k.  The  a  priori  model  probability  Pmm  t  is  therefore  given  by 

M 

Pmm,t  =  Pmm/k\ k-U  =  Y,Pi,m,tPm,,t(k-l\k-l)  (2) 

1=1 


The  model  transition  probabilities  are  a  priori  and  depend  on  the  target  model.  In  the  case  of 
track  initiation  on  a  sensor  with  Doppler  processing,  a  separate  target  model  is  defined  for  each 
velocity  ambiguity  and  there  is  no  transition  between  models.  However  if  the  models  represent 
different  target  manoeuvre  cases,  there  is  transition  between  models. 

2.1.2.  Clutter  Model 


The  clutter  is  modelled  as  uniform  in  one  or  more  measurement  dimensions  and  non-uniform 
in  the  remaining  dimensions.  At  least  one  dimension  is  assumed  to  have  a  non-uniform 
distribution.  It  is  assumed  that  the  measurements  are  due  to  more  than  one  pdf.  The 
surveillance  region  is  divided  in  the  uniform  dimensions  into  clutter  regions  each  of  which  has 
it's  own  non-uniform  pdf  in  the  remaining  dimensions.  There  are  C  different  non-uniform  pdfs 
and  hence  C  clutter  regions.  The  symbol  for  a  clutter  region  is  Oc  where  c  =  1, . .  ,C.  For  each 
measurement  y,  there  is  an  associated  clutter  region  index  c,.  The  number  of  measurements 
contained  in  a  clutter  region  is  Nc  and  its  area  (volume)  in  the  coordinates  for  the  uniform 
distribution  is  given  by  Ac.  The  sets  representing  all  clutter  regions,  measurements  and 

volumes  are  respectively  O  = 


The  coordinates  that  do  not  have  a  uniform  distribution,  such  as  signal-to-noise  (SNR),  are 
represented  by  a  continuous  clutter  distribution.  These  coordinates  are  identified  by  a  separate 
component  of  the  measurement  vector  so  that  y,  =  {z„  s,}.  Here  z,  contains  the  components  that 
have  a  uniform  clutter  distribution  in  the  clutter  region  with  index  c,  and  s,  contains  the 
components  with  a  nonuniform  distribution. 


Clutter  measurements  produced  by  the  sensor  are  assumed  to  be  independent  of  one  another. 
The  a  priori  probability  distribution  of  a  clutter  measurement  in  the  s  variables  given  it  is  in  a 
clutter  region  is  given  by  the  function  gc(s).  In  the  z  coordinates  the  distribution  is  uniform  so 

g  (s) 

that  the  pdf  of  clutter  measurements  is  given  by  g(y)  -  — — . 


3 


DSTO-TR-1184 


2.1.3.  Track  Data 


The  a  priori  track  data  consists  of  the  state  estimates  for  T  targets  where  each  is  represented  by 
M  models.  For  each  track  there  are  M  model  state  estimates,  covariances  and  model 
probabilities.  In  addition  the  track  data  includes  the  probability  of  target  visibility.  All  past 
data  up  to  time  k-1  are  given  by 


Z(fc- 1) 


{M-C, . 

{{*m/k^\k  V,-Pmt(k-l\k-l),  Pmmt(k-l\k-l)}^  ,  Pvt(k-l\k-l)y 


(3) 


’Dre  target  model’s  state  estimate  for  time  k  based  on  past  data  up  to  time  is  denoted  by 
*m,f  k  I  an<^  the  covariance  by  Pmt  ( k  \  j)  where  t  is  the  track  index  number  and  m  is  the  model 

index  number.  The  collection  of  M  estimates  is  called  track  t  and  the  target  associated  with  this 
track  is  designated  as  target  t.  It  is  assumed  that  all  targets  in  the  surveillance  area  have  an 
associated  track  representation.  The  state  estimate  of  target  t  is  denoted  by  xt(k\j)  and  the 

covariance  of  this  estimate  by  P  f  (k  \  j) .  These  are  obtained  from  the  sum  of  the  product  of  the 
model  state  and  probability. 


As  well  as  the  above  state  estimates  and  covariances,  the  track  data  base  also  has  the  target 
visibility  from  the  previous  time  k-1,  ie  Pvt(k-1  \  k-1 )  and  the  model  probability  Pmm  l(k-1  \  k-1)  for 

M 

the  M  target  models.  The  model  probabilities  are  constrained  by  Pmm  f  =  1  and  allowance 

m= 1 

is  made  for  the  model  to  change  with  the  time  index  k. 


The  model  transition  probability  matrix  gives  the  a  priori  probabilities  for  the  current  time  k. 
This  involves  combining  the  states  of  all  models  that  can  transition  to  this  model  from  time  k-1. 
After  some  maths,  the  a  priori  state  estimate  and  associated  covariance  conditioned  on  it  being 
model  m  is  obtained  by  the  following  equation. 


f  iVi 

nm,t  I  1 


^  M 

P™<*:|^1)  =  7^—  SI P,mJ ♦ 

nm,l 
1  M 

—  E  p,,4 W* I*-i>W*I*-dt  - 


(4) 


Pm  .  /.] 

m,t  1  1 


tm/k\k-Wm,W-l)r}Pmu(k-l\k-l) 


Where  xIM(.k\k-l)  -  Fm,  [fc-1)  and  .  Qm,. 
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For  the  special  case  of  track  initiation  where  velocity  ambiguities  exist  in  the  sensor 
measurements,  a  separate  target  model  is  allocated  to  each  velocity  ambiguity.  During 
initiation  the  ambiguity  does  not  change,  so  there  is  no  switching  between  models.  This  then 
makes  the  probability  transition  matrix  equal  to  the  identity  matrix  and  from  equation  (2), 
Pmmt(k\k-l)  =  Pmint(k-l\k-l).  So  equation  (4)  with  an  identity  matrix  for  the  probability 
transition  matrix  gives  the  standard  prediction  equations,  namely. 

T  (5) 

*V*I*-D  =  V^Cfc-llfc-DF^  -Q^ 


2.2.  Propagation  and  Receiving  Process 

The  model  for  the  propagation  and  receiving  process  in  Figure  1  employs  the  notion  of  target 
visibility,  ie  a  target  is  either  visible  or  not  visible  (invisible)  to  the  sensor.  Visible  targets  give 
sensor  measurements  that  satisfy  the  sensor's  target  and  measurement  models.  Conversely, 
Invisible  targets  give  no  measurements  that  satisfy  any  of  these  models. 

The  factors  that  determine  the  transition  between  the  visible  and  invisible  states  of  a  target  are 
independent  of  any  target  model.  Following  the  approach  in  [1],  the  a  priori  probability  that 
the  target  is  visible  is  given  by. 


Pvt  =  Pvt{k-\\k-l)Avt{k)  +  (l -Pv,  (k-l\k-l))A*(k) 
=  (a 7<fc)  -  A*(k))pvt(k-I\k-D  +  A f(jfc). 


(6) 


Where  Pvt  (.)  is  the  probability  of  target  visibility  for  target  t.  A  J(k)  is  the  a  priori  probability 

of  target  t  remaining  visible  from  time  k- 1  to  k,  and  A^(/c)  is  the  a  priori  probability  of  target  t 
changing  from  being  invisible  at  time  k-1  to  being  visible  at  time  k. 

For  notational  simplicity  Pvt  represents  the  a  priori  probability  Pvt(k  |  k-1).  When  denoting  the 
a  posteriori  probability,  the  time  index  is  included,  ie  Pv,(k\k). 

2.3.  Measurement  Process 

The  data  input  to  the  measurement  process  is  the  received  and  processed  signal  in  each 
resolution  cell.  The  signal  is  the  sum  of  the  returns  from  targets,  clutter  and  noise.  Clutter  and 
noise  are  collectively  called  clutter.  The  assumed  functions  of  the  measurement  process  are: 

(1)  Convert  the  receiving  system  output  into  measurements  that  contain  estimates  of 
parameters  such  as  range  and  bearing 
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(2)  Partition  the  receiving  system  output  into  the  C  clutter  regions  in  which  measurements 

from  clutter  are  represented  by  a  common  probability  density  function  (pdf). 
Determine  the  clutter  region  parameters  and  update  the  data  for  the  clutter  region  pdf 
models. 

The  details  of  the  above  steps  are  given  in[5]  [6]  and  [7]. 


The  measurements  output  from  step  (1)  are  represented  by  the  set  {y.}^,  where  y,  is  the  2th 

measurement  and  N  is  the  number  of  measurements  at  time  k.  Here  the  time  index  is  omitted 
for  clarity. 

Step  (2)  replaces  the  commonly  used  uniform  pdf  for  clutter  by  multiple  nonuniform  clutter 
regions.  It  determines  the  clutter  region  shape  0={0£}^=1/  and  the  clutter  region  volume 

A  =kVc=i  f°r  the  dimensions  where  the  clutter  is  uniformly  distributed.  For  each  measurement 

y>  R  determines  the  associated  clutter  region  index  c;.  It  also  calculates  Nc  the  number  of 
measurements  contained  in  the  clutter  region  Oc.  Using  the  measurements  yir  this  step  updates 
the  a  priori  distributions  for  the  clutter  regions. 

2.3.1.  Target  Probabilities 

The  target  probabilities  relate  to  when  the  target  is  visible  and  are  derived  from  target  signal 
fluctuations  and  the  conditional  tests  in  the  measurement  process.  The  measurements  are 
random  variables  and  are  independent  at  each  time  k  and  the  measurement  process  can  stop 
the  formation  of  a  target  measurement.  The  mechanism  that  causes  a  failure  to  form  a  target 
measurement  is  different  to  that  which  causes  a  target  to  be  invisible.  The  action  of  getting  a 
target  measurement  is  called  detection.  The  a  priori  probability  of  detection,  Pdmt,  is  the 
probability  of  obtaining  a  sensor  measurement  when  target  t  is  visible  and  obeys  the  222th  model. 

The  sensor  measurements  from  Visible  targets  are  related  to  the  target  state  by. 


y{k)  =  HmA®Xm/V  +  (7) 

Where  y(k)  is  the  measurement  vector,  H m  t(k)  is  the  measurement  matrix  for  model  m  of  track 
t  and  v{k)  is  the  measurement  noise  vector,  Gaussian  distributed  with  zero  mean  and  covariance 
R(k).  Here  H m  l(k)  and  R (k)  are  assumed  known  at  time  k.  The  model  index  for  H  is  included 
to  allow  for  cases  such  as  fixed  transponder  signals  in  Doppler  radars.  Here  the  measured 
Doppler  bears  no  relation  to  target  range-rate. 

The  pdf  for  target  measurements  is  defined  given  the  target  obeys  the  222th  model.  It  is  expressed 
with  respect  to  the  predicted  position  of  the  track.  The  distance  of  a  measurement  from  this 
position  is  normalised  by  the  innovation  covariance  for  the  222th  target  model  Sm  t(k),  where 

rfn/®  +  RW  •  The  equation  for  calculating  this  distance  is. 
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Where  y .  w ( is  the  vector  between  the  track  predicted  position  for  model  m,  ie  xf  m(k\k-l) ,  and 
the  z*  measurement,  ie 


Vim  =  yt,m,t  -  !*-!)• 


The  pdf  of  target  measurements  is  Gaussian  distributed  and  given  by. 


fm/y  i,m,) 


(2n)»/: 


(9) 


Where  p  is  number  of  degrees  of  freedom,  ie  the  number  of  components  in  the  measurement 
vector.  This  distribution  is  not  truncated  and  therefore  integrates  to  one  over  all  y. 

2.3.2.  Clutter  Probabilities 

The  clutter  density  is  correlated  with  time,  so  the  clutter  density  model  in  Figure  1  stores 
filtered  histograms  for  each  of  the  C  clutter  regions.  A  simple  low  pass  filter  for  each  clutter 
region  c  with  weight  wc(k)  is  applied  to  the  raw  histogram  data  from  the  measurements  ty,}^ . 

The  estimated  histogram  for  bin  b  is  given  by  [6]: 


=  hc  (k  - 1)  +  wc(k) 


h\k) 


ticb(k-l) 


GO) 


where  hb(k)  is  the  number  of  measurements  at  time  k  for  bin  b  and  clutter  region  c.  This 

histogram  data  is  used  to  compute  a  sample  cumulative  distribution.  This  distribution  is  then 
used  with  the  Kolomogorov-Smimov  test  statistic  to  estimate  the  clutter  density  gc(s)  as 
described  in  [6].  This  density  is  considered  a  priori  because  it  is  derived  from  filtered  past  data 
and  because  the  number  of  measurements  used  for  its  estimate  are  much  greater  than  the 
number  selected  by  any  track. 

2.4.  Measurement  Selection  Process 

The  M  models  for  each  track  t  select  a  fixed  number  of  I  measurements  from  the  set  {y,}  based 
on  their  distance  dtmi  away  from  the  predicted  model  position.  It  is  assumed  that  the 
probability  of  any  measurement  outside  of  the  I  nearest  being  caused  by  the  model  of  interest 
is  zero.  Using  a  fixed  number  of  neighbours  is  computationally  efficient  and  it  allows  the 
selection  region  to  expand  when  the  predicted  position  is  not  near  to  any  measurements  (as 
would  be  the  case  for  a  manoeuvring  target). 

Due  to  the  non-linear  relationship  of  the  MM-UJPDAF  computational  complexity  with  the 
number  of  target  models  and  tracks,  it  is  advantageous  to  divide  the  tracks  into  clusters.  We 
will  define  a  track  cluster  as  a  collection  of  measurements  and  tracks  that  are  assumed  to  form 
a  complete  set  for  probability  calculations.  Each  cluster  is  treated  individually  by  the  MM- 
UJPDAF.  Models  of  tracks  within  the  same  cluster  may  or  may  not  have  an  impact  on  the 
probability  calculations  for  the  track  model  of  interest,  but  it  is  assumed  that  all  tracks  outside 
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the  cluster  have  a  negligible  impact.  There  are  many  ways  that  these  track  clusters  may  be 
formed.  Ultimately,  the  size  of  the  clusters  is  governed  by  the  real-time  processing  constraints. 

The  process  of  forming  the  track  clusters  is  referred  to  as  clustering.  The  most  desirable 
clustering  method  is  to  group  all  tracks  whose  models  select  common  measurements.  This 
approach,  however,  leads  to  clusters  which  are  too  large.  Pruning  techniques  are  therefore 
applied  to  keep  the  cluster  size  manageable.  Because  of  the  pruning  performed,  clustering  must 
be  repeated  for  each  track  where  each  track  is  treated  as  a  reference  track  of  the  cluster. 

Clustering  for  this  filter  follows  the  basic  principle  that  only  tracks  which  are  likely  to  be  the 
true  source  of  measurements  selected  by  a  model  from  the  reference  track  are  added  to  that 
track's  cluster.  All  the  selected  measurements  for  each  model  of  each  track  in  the  cluster  are 
included  in  the  cluster.  This  approach  limits  the  cluster  size  to  a  more  manageable  level. 


Consider  the  example  illustrated  in  Figure  2.  Here,  tracks  are  denoted  by  crosses, 
measurements  with  circles  and  an  ellipse  has  been  drawn  to  show  which  measurements  are  the 
nearest  to  each  track.  In  this  example,  I  is  three.  There  are  two  models  for  each  track  (this  was 
chosen  for  simplicity  of  the  diagram).  Each  of  the  models  is  labelled  t(m )  indicating  model  m 
or  track  t.  When  we  form  the  cluster  for  track  1  we  include  track  2  but  not  track  3.  Clearly  this 
is  sub-optimal  since  track  3  will  impact  on  the  association  probabilities  for  track  2.  This 
approximation  is,  however,  necessary  in  order  to  calculate  the  probability  terms  and  it  is 
assumed  that  the  contribution  of  tracks  like  track  3  are  a  second  order  effect.  The 
approximation  is  particularly  necessary  in  high  target  density  conditions  where  many  tracks 
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select  the  same  measurements  and  when  the  JPDAF  is  incorporated  with  the  multiple-model 
initiation  filter. 

The  cluster  size  is  also  restricted  using  a  predefined  distance  threshold.  Tracks  are  only  added 
to  the  cluster  if  the  distance  between  the  track  and  a  peak  selected  by  the  reference  track  is  less 
than  this  predefined  threshold.  This  threshold  controls  the  cluster  size  and  is  a  design 
parameter.  The  distance  between  the  track  and  a  peak  is  the  same  normalised  distance  used 
for  measurement  association  and  so  can  be  viewed  as  a  threshold  measured  in  standard 
deviations  of  the  track  pdf. 

The  clustering  process  is  then  summarised  by  three  steps: 

•  Choose  a  reference  track 

•  Find  all  tracks  that  are  less  than  a  predefined  distance  away  from  any  of  the 
measurements  selected  by  the  reference  track 

•  Find  all  the  measurements  selected  by  any  track  in  the  cluster 

Once  the  cluster  has  been  formed,  it  is  common  [8]  to  form  a  Validation  Matrix  for  the  cluster. 
The  validation  matrix  shows  which  measurements  have  been  selected  by  each  track.  The 
matrix  can  be  used  to  determine  the  joint  events  [8].  The  approach  adopted  by  the  author  does 
not  make  use  of  a  validation  matrix,  but  instead  uses  track  lists.  These  lists  enable  the  event 
probabilities  to  be  written  in  a  factorised  form  and  hence  be  efficiently  calculated.  The 
Validation  Matrix  is  not  described  here  to  avoid  unnecessary  complication.  The  track  lists  are 
described  in  more  detail  in  the  following  section. 
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3.  Event  Space 

The  event  space  consists  of  all  possible  track  and  model  measurement  associations  for  a 
particular  track  cluster  at  time  k.  Define  IT  and  T  to  be  the  number  of  measurements  and  tracks 
in  the  cluster  respectively.  Let  y  refer  to  a  measurement  in  the  cluster  and  t  refer  to  a  track. 
Recall  that  the  number  of  models  for  each  track  is  M.  The  reference  track  is  track  1.  Each  model 
of  each  track  selects  I  measurements  as  possible  association  candidates.  These  I  are  the  closest 
measurements  to  the  a  priori  model  position  and  it  is  assumed  that  all  other  measurements 
have  a  negligible  probability  of  being  caused  by  this  particular  model. 

3.1.  Event  Space  Partitioning 

Define  the  joint  events  as  a  partition  of  the  event  space  such  that  within  that  partition  every 
measurement  in  the  cluster  has  a  unique  identified  source.  The  source  may  be  either  clutter 
(noise)  or  one  of  the  models  from  one  of  the  target  tracks  in  the  cluster.  No  more  than  one 
measurement  can  have  the  same  track  source  since  each  target  produces  at  most  one 

measurement.  The  joint  events  are  of  interest  because  it  is  possible  to  calculate  their 
probability. 

Let  the  symbol  0A  denote  a  particular  joint  event  which  defines  the  track  and  model  source  for 
each  of  the  IT  measurements.  The  subscript  A  is  the  joint  event  index.  Of  these  IT 
measurements,  nc  are  contained  in  the  clutter  region  c  and  nc  A  is  the  number  of  measurements 
assigned  to  targets  for  joint  event  A  in  this  clutter  region.  A  is  an  IT  element  list  whose 
elements  X.={m.,t^  are  indices  defining  the  model  and  track  that  are  the  source  of  the 

measurement  y}  under  the  particular  joint  event.  Let  X.  -  ( 0 , 0 )  represent  the  jth  measurement 

being  due  to  clutter  and  X.  ~{rn,t^  represent  the  fh  measurement  being  due  to  model  m  of  track 

t.  Each  track  element  f;  of  A  takes  a  value  in  the  range  [0,  T  ].  The  model  elements,  each  take 
a  value  in  the  range  [0,  M],  Since  each  target  may  cause  only  one  measurement,  no  two  track 
elements  of  A  can  have  the  same  non-zero  value.  In  fact,  the  range  of  possible  values  for  Xf 
is  also  limited  by  the  measurement  selection  process.  Since  each  model  may  only  consider 
associations  with  the  I  nearest  measurements,  Ay  can  only  take  the  value  ( m,t )  if  measurement 
y> 1S  one  of  the  1  nearest  measurements  selected  by  model  m  of  track  t.  To  quantify  which 
models  from  each  track  select  each  measurement  two  further  lists  are  defined. 

Define  Tj  as  the  number  of  tracks  selecting  measurement  y-y  To  describe  the  possible  values  of 
Xj,  define  the  track  list  t.  =  jt.(f)  whose  elements  are  the  tracks  which  have  selected 

measurement  yy  Let  the  first  element  t;(0)  be  zero  (denoting  that  clutter  is  a  possible  source  of 
the  measurement)  and  element  numbers  1..7}  take  values  in  [1,  T]  corresponding  to  the  tracks 
which  have  selected  yf 

The  symbol  Mj  t  defines  the  number  of  models  from  track  t  that  have  selected  measurement  yr 
Mj ,t  is  zero  for  all  tracks  not  in  the  list  t ).  Define  the  model  list  p  =  jp.  (m)  \M'-‘  whose 

elements  are  the  models  from  track  t  that  have  selected  measurement  yr  For  clutter,  let 
^/,o  ~  {0}  with  M.q  =  1  for  all  measurements.  These  lists  are  required  when  we  wish  to 
enumerate  the  joint  events. 
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3.2.  Marginal  events 

Based  on  the  solution  of  the  MM-UPDAF  filter  [1],  the  MM-UJPDAF  solution  requires  the 
probability  of  each  of  the  7  +  2  events  for  each  model  of  the  reference  track.  These  events  are 
composed  of  the  joint  events  described  above  and  are  referred  to  here  as  marginal  events.  This 
name  is  used  because  these  events  are  a  union  over  the  joint  events  analogous  to  forming  the 
marginal  pdf  for  one  variable  of  a  multi-variable  joint  pdf. 

Denote  the  z-th  marginal  event  for  the  227th  model  of  the  reference  track  as  0,m.  Let  the  pointer 
k (z',222)  identify  the  2th  closest  measurement  to  the  222th  model  of  the  reference  track.  That  is,  yK(j  m) 
is  the  2th  closest  measurement  to  the  222th  model  of  the  reference  track.  For  notational  brevity,  the 
dependence  of  k  on  the  model  m  and  the  index  2  is  dropped  such  that  k  =  k(z,22z)  and  yK  s  yK(i  m). 

The  7+2  marginal  events  for  each  model  of  a  reference  track  are: 

©_!  m  the  model  m  is  the  correct  model  and  the  target  is  not  visible 
©0  m  the  model  m  is  the  correct  model  and  the  target  is  visible  but  none  of  the  7 
measurements  selected  by  it  are  caused  by  the  track 
©;  m  the  model  m  is  the  correct  model  and  the  target  is  visible  and  has  caused  the  ilh 
closest  measurement  to  it  (that  is  measurement  yK).  Here  z  must  lie  in  the  range 
[1,7] 

To  derive  an  estimate  of  the  track  for  a  particular  model  m,  the  reference  track  model  event  Qm 

2 

is  defined.  This  is  given  by  Qm  -  U  0.  . 

2  =  -1  t,m 


3.3.  Relationship  between  Joint  and  Marginal  events 

For  the  marginal  event  0,  m ,  model  m  from  the  reference  track  is  the  source  of  measurement  yK 
when  2  lies  in  the  range  [1, 7  ].  This  event  is  the  union  over  all  joint  events  where  model  m  from 
the  reference  track  is  the  source  of  measurement  yK.  That  is,  all  joint  events  where  XK  =  (222,1). 


The  union  relationship  is  mathematically  expressed  as: 


©.. 


U 

A  IV  (m,l) 


0A 


(11) 


where. 


U 

A  |XK=  (m,l) 


uu, ...  uu  .1x1 


K-l 


XX  eA 


IT 


with 


t/f; 

U 

u 

f  -v® 

(12) 


(13) 


and 
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K'*i  ) 

fac-A-i) 


AH 


(mjr'  ^rr) 


(14) 


The  modified  track  lists  are  defined  by 

Vfoiu^nnvy)  as) 

and 

f .  =  card  (t.)  -  1  (16) 

The  modified  track  list  x.  is  simply  the  track  list  for  measurement  y;  without  the  tracks  that 

have  already  been  associated  with  a  previous  measurement.  Using  the  modified  track  lists 
ensures  that  no  track  is  associated  with  more  than  one  measurement  for  any  particular  joint 
event. 


Note  that  the  union  denoted  UU.  is  simply  an  enumeration  of  all  the  track  models  that 

measurement  y}  can  associate  with,  given  the  associations  already  made  by  measurements 
1  .  .  j- 1.  The  order  of  the  unions  cannot  be  interchanged  because  of  this  dependance  on  the 
associations  with  previous  measurements.  This  is  reflected  in  the  recursive  nature  of  the 
measurement  modified  track  list  f .. 

i 

Since  the  joint  events  are  mutually  exclusive,  it  is  apparent  that  the  posterior  probability  of  a 
marginal  event  is  simply  the  sum  of  the  posterior  probability  of  the  appropriate  joint  events. 

The  events  corresponding  to  i  =  0  and  i  =  -1  must  be  treated  slightly  differently.  The  reason  for 
this  is  that  the  joint  events  cannot  differentiate  between  the  two.  The  joint  event  index 
determines  that  the  reference  track  has  caused  no  measurement  but  it  does  not  describe  which 
model  is  correct  or  whether  the  target  is  visible.  The  union  of  all  joint  events  where  the 
reference  track  causes  no  measurement  is  then  equivalent  to  the  union  over  all  models  of  the 
events  &0  m  and  @_x  m.  The  expression  for  this  union  is  given  below: 


M 

Ue, 

m=l 


U 

A  \l^t  *  (m,l)  \/m,j 


(17) 


where. 
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U  0A  =  W,  ...  LU;t  0a 

A  \Xj  Mm,j 


(18) 


and 

("V  h  ) 

A=  ) 

(m/r  trr) 

The  modified  indices  are  the  same  form  as  previously  presented  in  equation  (15)  and  the  unions 
UU.  are  as  presented  in  equation  (13) 

Once  again,  the  posterior  probability  of  the  union  of  Q.j  m  and  @0  m  is  calculated  by  summing  the 
posterior  probabilities  of  the  relevant  joint  events.  To  calculate  their  individual  probabilities, 
the  conditional  probabilities  Pr{0  a  J©.,  ^©^IL-ll©  1MU©0M}  anc* 

Pr  {0o,m  1 0-i,i U  0o, i U ■ "  U  0-i ,mU  0o,m^  are  squired. 

3.4.  The  Joint  Event  Tree 

One  way  to  better  understand  the  operation  of  equations  (12)  and  (18)  is  to  use  a  graphical  aid 
to  visualise  the  event  enumeration.  It  is  appropriate  to  use  a  tree  structure  as  this  graphical  aid. 
The  formation  of  this  joint  event  tree  is  illustrated  by  way  of  example. 

Consider  a  cluster  that  contains  three  tracks  and  four  measurements.  For  the  moment,  let  each 
track  have  only  one  model.  The  measurement  selections  of  the  tracks  are  described  by  the  track 
lists  below. 

Tj  =  {0, 1,  2}  with  T1  =  2 

t2  =  {0, 1,  2, 3}  with  T2  =  3 

t3  =  {0, 2,  3}  with  T3  =  2 

t4  =  {0, 1 ,3,}  with  T4  =  2 

In  this  example,  1  =  3,  IT  =  4  and  T  =  3 

Let  us  now  consider  the  case  where  track  1  causes  no  measurement,  ie  equation  (17)  applied 
to  this  cluster.  Starting  with  the  first  list,  =  {0,1,2},  equation  (15)  is  used  to  determine  x1  as 

follows. 

x,  =  {0}u{  {0,1,2}  nul} 

=  {0}  U  { {0,1,2}  n  {0,2,3} } 

=  {0,2} 
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Since  the  first  measurement  is  not  selected  by  track  3,  only  clutter  and  track  2  are  considered. 
These  values  are  then  taken  in  turn  to  generate  the  modified  track  lists  ?2 .  Using  equation  (15) 

as  above,  the  lists  and  the  associated  value  for  t  are: 

|  {0,2,3}  when  t1  =  TjfO)  =  0 
|  {0,3}  when  ^  =  t^l)  =  2 

A  simplified  way  to  undertake  these  steps  and  determine  the  track  lists  is  to  use  the  joint  event 
tree  shown  in  Figure  3.  Each  node  on  this  tree  gives  the  track  associated  with  a  measurement 
and  each  layer  corresponds  to  a  particular  measurement.  The  leaf  nodes  of  the  tree  contain  the 
source  of  the  measurement  yIT  and  the  path  through  the  tree  to  reach  the  leaf  node  indicates  the 
source  for  each  of  the  other  measurements.  The  children  at  each  node  represent  the  list  of 
possible  track  associations,  x.. 


The  example  above  shows  that  the  branches  available  for  selecting  the  value  of  t2  depend  on 
what  value  t1  has  taken.  This  is  also  highlighted  by  the  modified  track  lists  preceding  Figure 
3.  From  the  joint  event  tree,  the  left  most  path  gives  A  =  {0, 0, 0, 0} ,  ie  all  measurements  are  from 
clutter.  Likewise  the  rightmost  path  gives  A  =  {2, 3, 0,0} . 

Adding  more  models  to  the  tracks  does  not  change  the  tree  structure.  When  there  is  more  than 
one  model,  each  node  on  the  tree  represents  a  track  choice.  After  the  track  is  chosen,  a  model 
from  that  track  is  chosen  using  ji/  (  to  determine  which  models  have  selected  that  measurement. 
The  models  add  breadth  to  the  tree,  however,  since  each  track  may  only  have  one  true  model, 
the  lower  branches  of  the  tree  are  the  same  for  each  model 

The  joint  event  tree  is  a  useful  tool  not  only  because  it  presents  a  formalised  exhaustive 
enumeration  scheme  but  also  because  of  the  factorisation  it  represents.  Using  the  tree  type 
approach  to  enumerating  the  events  is  much  like  the  divide  and  conquer  strategy  used  in  such 
algorithms  as  fast  fourier  transform  and  quicksort  [10].  This  factorisation  is  important  in  the 
calculation  of  the  marginal  event  probabilities  as  described  in  the  following  section. 
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4.  Target  State  Estimate 

The  posterior  state  estimate  for  the  JPDA  filter  is  the  same  as  that  used  for  the  single  target 
multiple  model  PDA  filter  in  ref  with  the  exception  of  the  event  probability  calculations. 
However,  the  state  estimate  is  repeated  here  for  completeness  and  clarity.  To  simplify  the 
equations,  the  track  index  t  is  omitted. 

From  Section  2,  the  data  provided  by  the  sensor  with  each  multiple  model  track  consists  of  the 
selected  measurements  Y(k),  the  regions  containing  the  I  nearest  measurement  SI,  the  clutter 
region  parameters  at  time  k,  ie  n,  N,  O  and  the  past  track  data  Z(k- 1).  For  notational  clarity  the 
time  dependence  of  the  sensor  data  is  not  shown.  From  the  defined  target  and  sensor  model 
in  Section  2  and  the  marginal  events  defined  in  Section  3.2,  the  optimal  estimate  of  the  target 
state  x(k)  is, 

x(k\k)  =  E  {x(k)  |  Y,  SI,  n,  N,  O,  Z] 

m  i  (20) 

E  E  vl„rnxm(k)\elm,Y,si,n,N,o,z}. 

m-1  i--l 

Where  Pim  =  Pr{0.jY  ,SI,n,N,0,Z} 

0.  =  is  the  marginal  event  as  defined  in  Section  3.2 

The  evaluation  of  E{xm(k)  |  0.  m,  Y,  SI,n,N,  O,  Z }  for  i  >  0  is  the  Kalman  filter  solution  for  updating 
the  predicted  state,  xm(k  |  /c-1) ,  of  a  target  represented  by  model  m,  with  measurement  yi .  Thus, 
equation  (20)  may  be  written  as  the  following  summation. 

M  I 

#*|»  =  E  E  <21> 

m= 1  i=-l 

x  (k\k-l)  +  W  y.  i  =  l,..,/ 

xm{k\k-l)  i  =  -1,0 

-  Q>Tsm' 

ie  the  Kalman  gain  matrix  for  the  rtf'  model. 

From  equation  (21),  each  model  is  similar  to  a  PDA  filter.  These  remain  separate  with  the 
contribution  of  a  model  to  the  state  estimate  tending  to  zero  as  its  probability  approaches  zero. 
This  filter  therefore  combines  the  approaches  taken  in  PDA  and  track  split. 

The  selection  process  requires  the  estimate  of  each  model’s  state  given  it  is  valid.  This  is 
obtained  by  writing  equation  (20)  in  the  following  form. 
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Where 


x(k\k) 


*<*!*)  =  £  PJm(k\k). 


Pr{Qm| Y,SI,n,N,0,Z}  =  £  p. 

i=-l  ' 

E{x  (fc)|Q  ,Y,S/,rc,N,0,Z}. 


From  equations  (20)  and  (22),  the  estimate  of  each  model's  state  is  then  given  by, 

*«<* l*>  =  *m<k|*-l)  +  Wwy~  (fc|fc). 


Where 


Now  the  covariance  of  £(k\k)  is  a  weighted  sum  of  Gaussians  where  each  is  derived  from  the 
PDAF  approximation  for  each  model.  The  event  probabilities  for  each  model  require  the 
covariance  for  £Jk\k)  which  is  denoted  by  Pm(k\k).  This  is  derived  in  Appendix  A  and  is  given 
by  the  expression. 

VJklk)  =  £i=P. lm(*|fc-l,  *  ♦  (  1  a  ft|t) 

/  x  (24) 

1  (3  \ 

+  W  Y'  — —  y.  vT  ~  v  v1  WT 

,  M  Pm 


Where  Sm(^r  |  j)  is  the  Kalman  filter  covariance  for  model  m  at  time  k  given  data 

to  time  /. 

a  is  the  covariance  for  model  m  when  the  target  is  invisible  at  time 
k  given  data  to  time  j. 

The  display  of  track  position  to  an  operator  requires  the  state  estimate  £(k  |  k)  while  multi-sensor 
track  fusion  uses  the  covariance  P(*|*).  This  is  obtained  from  the  model  covariance 
Pm(£|A:)using  the  following  expression. 

M 

PttlW  -  E  PJP.ttlW  -  «„  -»«.  -  »T]  (25) 
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The  derivation  of  the  event  probabilities  P;  m  required  for  equations  (20)  and  (24)  is  given  in 
Appendix  A.  These  probabilities  simplify  to, 

b._ 


p.  - 

M  / 

E  E  b, 

M/-1 


(26) 


Where 


Km  =  E  MVyt,) 

ty^O) 


E  t(yK+A+ ,) 

Cl=W0> 


-  E  ^k-A-i) 


f  _  . 

\ 

- 

t,t(T(7.)  C 

( N  -n 

... 

E  'KyJT,f/T)  n 

c  c,  A 

fl  -fl 

... 

^  t C- 1 

\  c  c'K) 

/ 

b0,m 


Pm .  PvAl-Pd,  Ps ,  ) 

l,m  1  '  l,m  y' 

PP 


/  _  . 

C 

f  N  -n  .) 

E  MyIT,tlT)  n 

c  c,  A 

M  -ft 

^  l  ff~ ^/7*(0)  C-l 

^  c  c,A/ 

Pm ,  (1  -  PiO  W 

b -i,m  =  -  -  ,mi- -  E 

Pf,  *!=*!©) 


f  /  . 

\  j 

T/t(Tj7-)  C 

f  N  -n  .'j 

E  ^(y/r^/r)  IT 

c  c,A 

ft  ~Yl 

...1 

^  C  =  1 

\  c  c,K) 

)  ) 

in  the  above  equations. 


M  r  i 

Pf.  =  rPm.  <  1  -Ps.  Pd.  Pv.\ 

t  ifm  I  *  i,m  i  \ 

m=l 


I'mW  = 


Pdm,tPmm,tPv>fm,tW 

gc(y)ipit 


me[l,M],  fe[l,T] 

ra=0,  f=0 


s  E  'UyI*9 

m=Py(0) 

The  b  equations  above  are  somewhat  complex  and  are  more  easily  understood  by  revisiting  the 
joint  event  tree.  Consider  the  calculation  of  b0  m  for  the  cluster  example  given  in  section  3.4.  The 
joint  event  tree  for  this  cluster  is  shown  in  Figure  3.  To  calculate  the  event  probability,  we  use 
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the  joint  event  tree  and  place  values  on  each  node. 
Pmm,iPv,  «-pdmAPsmA) 


The  root  node  takes  the 


value 


Ph 


and  the  leaf  nodes  take  the  values 

c=l  v 


N  -n  . 

c  c,  A 
n~nc,A) 


The 


intermediate  nodes  take  a  value  of  i |/(y  .,h) .  This  new  tree  is  shown  in  Figure  4.  The  event 

probability  is  then  calculated  by  summing  the  leaf  nodes  and  multiplying  the  sum  by  the  value 
of  the  parent  node.  The  parent  node  value  is  replaced  by  this  product  and  the  leaf  nodes  are 
discarded.  This  process  is  repeated  until  all  the  nodes  of  the  tree  have  been  removed  and  the 
resulting  value  is  the  un-scaled  event  probability  b0  m 


In  Figure  4  only  the  left  branch  of  the  tree  is  shown  to  simplify  the  diagram.  The  clutter 
associations  have  also  been  filled  with  the  value  i|/(y/r0)  =1 .  The  number  of  clutter  zones  is 
assumed  to  be  one. 


..l  i  1 . 


[  My**)) 


!  1 


\*iyai\ 


[  Wy^)) 


l 

. .»* 


x  X  ,/  X  ...1 . /  i 

""  .  * .  * . **  * . * 


Figure  4.  Joint  event  tree  with  node  values 


X) 


i#y*2)} 


m 


The  a  priori  visibility  Pv  is  given  by  equation  (6).  From  the  event  probabilities  p.  m  the  a  posteriori 
target  visibility  is  given  by. 


Pv(k\k)  =  1  - 


M 


E  P-i,, 


(27) 


This  visibility  is  then  used  in  equation  (6)  as  the  past  data  in  the  next  time  step  with  the  a  priori 
transition  probabilities. 
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5.  Filter  Performance 

The  MM-UJPDAF  filter  has  been  implemented  to  operate  on  the  Jindalee  Over-the-Horizon 
Radar  (1RSU).  The  implementation  of  this  filter  differs  from  the  standard  and  most  intuitive 
method.  It  factorises  the  event  probability  equations  to  give  a  considerable  computation 
advantage.  The  advantages  of  the  factorised  approach  are  now  considered.  This  is  then 
followed  by  examples  of  performance  on  OTHR  data. 

5.1.  Computation  Load 

The  brute  force  method  of  implementing  the  JPDAF  is  to  enumerate  all  possible  joint  events 
and  calculate  their  event  probabilities.  When  the  event  probability  is  calculated,  the  algorithm 
determines  which  marginal  events  this  joint  event  contributes  to  and  accumulates  their  values. 
As  has  been  described  in  section  4  this  is  not  the  way  that  the  JPDA  algorithm  for  1RSU  has 
been  implemented.  Instead,  by  factorisation,  a  closed  expression  for  each  marginal  probability 
has  been  found  and  the  probability  is  calculated  directly.  This  leads  to  a  reduction  in  the 
computation  requirement  and  the  functional  complexity. 

There  are  a  number  of  existing  'fast  JPDAF'  algorithms  which  have  been  designed  to  reduce  the 
computational  requirements  of  the  filter.  These  algorithms  generally  take  a  similar  form  to  the 
brute  force  approach  but  intelligently  use  previous  event  calculations  to  reduce  the  effort 
required  to  calculate  the  probability  of  similar  joint  events.  One  such  algorithm  is  the  LSA  of 
[9].  This  algorithm  takes  advantage  of  factorisation  to  compute  the  joint  event  probabilities. 
The  LSA  is  a  more  efficient  method  than  the  MM-UJPDAF  for  calculating  the  marginal 
probabilities  in  the  absence  of  multiple  models.  However,  the  complexity  of  the  LSA  rises 
sharply  when  multiple  models  are  introduced. 

The  reason  why  the  LSA  is  more  efficient  than  the  presented  form  of  the  MM-UJPDAF  is  that 
the  MM-UJPDAF  forms  a  separate  cluster  for  each  track  and  calculates  the  joint  events  for  the 
cluster.  When  several  tracks  form  the  same  cluster,  they  will  share  common  joint  events, 
however  the  MM-JPDAF  does  not  take  this  into  account.  Remember  that  the  cluster  rules  for 
the  MM-UJPDAF  limit  the  cluster  for  the  case  of  a  reference  track.  When  these  constraints  are 
relaxed,  larger  cluster  sizes  are  likely  with  shared  common  joint  events.  It  is  difficult  to 
compare  the  event  calculation  because  of  the  difference  in  clustering.  However,  we  examine 
the  LSA  load  with  tracks  and  models. 

A  bound  on  the  number  of  multiplication  operations  required  for  the  core  of  the  LSA  method 
is  given  by  [9]  as: 


N™,,<  ECrL7L-’(/^-2)  (28) 

2 

assuming  that  each  track  selects  I  measurements. 

When  each  track  consists  of  M  models,  this  expression  becomes: 
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^ mult  <  ]C  (IxM)1  ’(/xM+L-2)  (29) 

L=2 

The  above  bound  is  loose  but  it  is  significant  to  notice  that  the  computational  complexity  is 
polynomial  of  order  IxM  ‘  \  For  a  particular  cluster  size  (fixed  I  and  T)  the  number  of 
computations  required  grows  at  a  rate  determined  by  an  order  T-l  polynomial  in  M. 

The  MM-UJPDAF  approach  requires  a  separate  tree  calculation  for  each  model,  however  each 
tree  is  independent  of  the  number  of  models  since  the  tree  nodes  contain  values  which  have 
been  summed  over  all  models  for  the  corresponding  association.  Also,  models  from  a  track  are 
likely  to  share  measurements.  Since  the  root  of  the  tree  is  the  only  point  dependent  on  the 
model,  tree  calculations  for  common  measurements  can  be  shared.  The  result  is  that  the  MM- 
UJPDAF  algorithm  grows  less  than  linearly  with  the  number  of  models. 

It  is  infeasible  to  compare  the  measured  computation  performance  of  the  MM-UJPDAF  with 

existing  approaches  since  considerable  effort  would  be  required  to  implement  these  alternative 
algorithms. 

5.2.  Performance  examples 

If  the  number  of  tracks  in  the  cluster  is  only  one,  then  the  marginal  event  probabilities  simplify 
to  the  single  target  filter  equations  as  presented  in  [1]  (as  one  would  expect).  Having  observed 
this,  it  is  apparent  that  the  performance  of  the  JPDAF  is  the  same  as  the  single  target  PDAF  for 
cases  where  the  target  tracks  are  well  separated  and  the  clustering  process  does  not  add 
neighbouring  tracks  into  a  cluster.  It  is  useful  then  to  only  consider  the  performance  of  the 
filter  in  circumstances  where  neighbouring  tracks  do  interfere  with  each  reference  track. 

5.2.1.  Multi-path  Performance 

In  OTHR,  the  most  common  source  for  closely  spaced  target  returns  is  the  multi-path 
propagation  medium.  Under  most  conditions,  the  ionosphere  provides  at  least  two  possible 
propagation  paths  between  the  radar  and  the  target.  In  the  simplest  model,  the  ionosphere  can 
be  considered  as  multiple  reflecting  planes.  The  two  main  propagation  layers  are  called  the  E 

and  F  layers.  With  these  two  layers  present,  there  exist  four  propagation  paths  as  shown  in 
Figure  5. 
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Figure  5.  Multi-path  propagation  in  OTHR 


The  radar  measures  range  using  the  measured  time  delay  so  that  different  paths  appear 
separated  in  range  due  to  the  different  path  length.  The  closest  paths  are  the  E-F  and  F-E  paths, 
called  the  mixed  mode  paths.  The  paths  both  give  the  same  time  delay  (that  is  group  or  slant 
range)  but  because  they  are  received  at  different  elevation  angles  the  coning  effect  of  a  linear 
array  causes  them  to  be  measured  at  different  azimuths.  The  elevation  of  the  transmitted  signal 
also  causes  slightly  different  radial  velocity  measurements.  Now,  since  the  radar  does  not  have 
a  measurement  of  tangential  speed,  the  tracks  must  be  initialised  with  fairly  large  azimuthal 
covariances.  The  nature  of  the  ionosphere  also  tends  to  make  the  azimuthal  measurement  more 
noisy  than  the  range  and  radial  velocity  measurements.  These  factors  conspire  to  produce  a 
tracking  environment  where  it  is  common  to  find  parallel  tracks  with  the  same  range  and  speed 
but  separated  in  azimuth.  Due  to  various  environmental  factors,  these  propagation  paths  tend 
to  fade  and  return  with  the  result  that  a  single  target  tracker  may  not  adequately  track  these 
paths. 

The  JPDA  has  been  demonstrated  to  improve  tracking  under  this  condition.  An  example  of  this 
improvement  is  shown  in  Figure  6.  The  single  target  PDAF  performance  is  shown  in  the  left 
plot  and  the  multi-target  JPDAF  performance  is  shown  on  the  right.  The  vertical  axis  is  range 
and  the  horizontal  axis  is  azimuth.  Here  the  single  target  tracker  fails  to  correctly  track  the  two 
mixed  mode  target  returns.  The  tracks  formed  by  the  single  target  filter  cross  over  when  the 
propagation  modes  fade  and  the  tracks  are  lost.  In  contrast,  the  JPDAF  holds  the  two  tracks 
along  their  separate  traces. 
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Figure  6.  Multi-path  tracking  performance 


5.2.2.  True  Multi-Target  Performance 

The  filter  also  improves  the  tracking  performance  for  more  traditional  multi-target  tracking 
scenarios.  Figure  7  below  shows  the  improved  performance  for  the  case  of  two  crossing  target 
tracks.  In  this  example,  the  two  tracks  have  similar  radial  velocities,  but  their  tangential  speeds 
are  different.  When  the  two  tracks  cross,  the  single  target  filter  drops  one  track  and  the  other 
track  switches  to  the  trajectory  of  the  dropped  track.  The  multi-target  filter  correctly  follows 
the  crossing  trajectories.  As  in  Figure  6  the  single  target  filter  is  shown  on  the  left  and  the 
multi-target  filter  on  the  right.  The  tracks  are  plotted  as  range  against  azimuth. 
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The  third  example  of  the  filter  performance  is  shown  in  Figure  8.  This  is  a  much  more 
congested  target  situation.  In  this  case  there  are  three  main  targets  each  being  propagated  over 
number  of  modes.  These  targets  are  moving  at  slightly  different  speeds,  but  they  are  very  close 
geographically.  In  addition  to  these  targets,  there  are  two  slower  targets  at  close  range.  These 
two  targets  are  actually  the  tracks  illustrated  in  Figure  7.  Due  to  the  dense  traffic,  it  is 
uninformative  to  view  the  tracks  on  a  range-azimuth  plot  as  in  Figure  6  and  Figure  7.  Instead, 
the  tracks  are  presented  in  Figure  8  as  range  against  time.  Each  of  the  propagation  modes  can 
be  observed  at  a  different  range  and  for  each  mode  there  should  be  three  tracks  -  one  for  each 
of  the  targets.  As  can  be  seen,  the  JPDAF  improves  the  performance  on  this  data  considerably 
and  separates  the  individual  targets  well.  The  single  horizontal  trace  at  lower  range  is  a 
synthetic  calibration  signal  injected  at  the  radar  receiver. 


Time 

Figure  8.  Single  Target  and  Multi-Target  Performance  for  High  Density  Targets 
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6.  Conclusion 

This  paper  has  developed  the  theory  for  the  multiple  target  extension  to  a  multiple  model  filter 
that  was  originally  developed  for  track  initiation  with  an  OTHR.  Including  multiple  targets 
with  multiple  models  greatly  complicates  the  evaluation  of  the  event  probabilities.  To  simplify 
this  step,  the  event  space  was  developed  with  the  aid  of  a  new  joint  event  tree.  The  event  tree 
allows  the  efficient  real-time  evaluation  of  the  event  probabilities  by  factorizing  the  event  space. 
This  new  approach  is  more  efficient  than  current  methods,  particularly  when  dealing  with 
multiple  models.  This  system  has  demonstrated  real-time  operation  when  the  track  cluster  size 
reaches  10  with  over  4  models. 


pie  multiple  target  extension  incorporates  all  the  previous  features  in  the  multiple  model  filter, 
ie  target  visibility,  multiple  non-uniform  clutter  zones  and  selection  of  measurements  based  on 
nearest  neighbours.  The  new  filter  was  tested  on  Jindalee  OTHR  recorded  data.  It  was 
implemented  for  the  tracking  stage  following  track  initiation.  Once  a  track  became  established, 
a  JPDA  filter  was  used.  This  tracking  system  demonstrated  enhanced  performance  when 
targets  become  closely  spaced.  With  distant  targets,  the  same  performance  is  obtained  as  for 
the  original  filter  without  multiple  targets.  So  this  new  filter  extends  the  performance  envelope 
of  the  previous  filter  and  differences  are  only  noted  under  conditions  of  closely  spaced  targets. 


Future  work  is  planned  to  add  the  multiple  model  filter  extension  to  the  JDPA  filter  for  tracking 
manoeuvring  target  during  established  target  tracking.  This  will  be  considered  with  other 
possible  enhancements. 


l 
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Appendix  A  Derivation  of  Event  Probabilities 


To  calculate  the  track  posterior  state  estimate  we  require  the  posterior  probabilities  of  the 
marginal  events.  These  probabilities  are  calculated  from  the  sum  of  joint  event  posterior 
probabilities  (refer  to  Section  3).  To  derive  an  expression  for  the  marginal  events,  we  begin  by 
seeking  an  expression  for  the  posterior  probability  of  an  arbitrary  joint  event.  ' 

Firstly  let  us  revisit  the  definitions  of  the  cluster  and  the  clutter.  There  are  IT  measurements 
and  T  tracks  in  the  cluster.  Each  track  has  M  models.  Let  IT  denote  the  whole  cluster.  The 
clutter  zone  associated  with  measurement  yf  is  c;.  O  is  the  set  of  clutter  zones  each  of  which 
contains  Nc  measurements  respectively.  Of  these,  nc  are  found  in  the  track  cluster  and  ac  of 
these  are  associated  with  tracks  under  the  event  0A. 


The  posterior  probability  of  the  a  joint  event  can  be  expanded  using  Bayes'  Rule  to  give: 


p(eA|i7,o,z)  = 


p(qa,u\z,o) 

P{u\z,o) 


p(dA,u\z,o) 

£p(eA,u|z,o) 

0A 


(Al) 


The  denominator  of  the  above  expression  can  be  considered  as  a  normalising  factor  and  is  the 
sum  of  the  numerator  over  the  whole  event  space.  We  need  therefore  only  consider  the 
numerator.  Recognise  that  the  cluster  consists  of  the  tracks,  the  measurements,  the  number  of 
each  of  these  and  the  track  selection  boundaries.  Since  the  track  positions  and  the  number  of 
tracks  is  fully  determined  by  the  past  data  these  features  are  not  probabilistic.  Thus  U  is 
expanded  into  the  measurements  Y,  the  number  of  measurements  n  and  the  cluster  boundary 
The  numerator  is  a  joint  probability  that  can  be  expanded  as  a  product  of  conditional 
probabilities  (again  using  Bayes'  Rule). 


p(eA,Lf|Z,o)  -  p[dA,Y,n,S\Z, o) 

=  P(Y\n,QA,S,Z,o)  p(n\0A,S,Z,O )  p(eA|S,Z,o)  P{S\Z,0) 


(A2) 


The  first  term  is  the  probability  of  the  measurements  given  the  number  of  measurements  and 
the  source  of  each  measurement  as  defined  by  the  joint  event.  Since  the  measurements  are 
independent,  this  is  simply  the  product  over  all  measurements  of  the  measurement  density  for 
each  measurement.  Thus,  the  first  term  in  equation  (A2)  is  given  by: 
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p(Y\n,s,6A,z,o)  =  n  fM  n  gjy? 

\eA\tfO  '  A..eA|f  =0  1 


fx  (?//)  IT 

n  iu.<y,> 

VAl tf0  8c,\/i'  H  ' 


(A3) 


where 


gc  (y^.)  is  the  probability  density  function  for  measurement  yj  being  from  clutter 
given  it  is  within  C  fl  S, 


ie 


,  with  Pc 

ci 


f  sM  dV 

yeS 


fx. 


(y;)  is  the  probability  density  function  for  a  measurement  caused  by  model  m, 


of  target  t,  visible,  detected  and  within  the  region  containing  the  I  nearest 
measurements  SI, 


A  (V 

ie  — - ,  with  Ps. 

Ps ,  tm 


f  f,Jy^dy 

yeS 


The  second  term  is  equation  (A2)  is  the  probability  of  observing  nc-ac  measurements  from 
clutter  zone  c  in  the  selection  area  given  the  total  measurement  density  and  the  size  of  the 
selection  area  defined  by  Nc,  Ac  and  ac.  This  is  given  by  the  binomial  distribution  and  is  written: 


'N  -n  an 

c,  A 


j>(nieA,s ,z,o)  -  n  „c_V'A  Pc'r’“ {'-pcc)"rn‘ 

c= 1  \  ,lc  c, A/ 


(A4) 


The  third  term  in  equation  (A2)  is  the  prior  probability  of  the  event.  This  is  the  prior 
probability  that  each  track  causes  the  output  described  by  the  event  index  A.  Since  the  tracks 
are  independent,  this  is  the  product  over  all  tracks  of  the  prior  probability  for  that  track. 


For  the  tracks  which  cause  measurements  under  A  we  require  the  prior  probability  that  Xl  m  is 
the  correct  model  for  track  Xl  t  and  that  this  track  is  visible,  detected  and  has  formed  the  ith 
measurement  inside  S.  This  probability  is  given  by  the  expression: 


where. 


PSm.t  Pdm,t  Pmrn,t  Pvt 


I 


ps„,  =  *K,i 


Pv,  =  p(v,  |z) 
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For  the  tracks  that  cause  no  measurement,  this  is  the  sum  over  all  models  of  the  probability  that 
model  m  is  the  correct  model  for  track  t  but  that  the  track  is  either  invisible,  not  detected  or  did 
not  cause  any  of  the  I  nearest  measurements.  This  quantity  is  given  by: 

M 

£  Pm,n,t{  ^~Ps,J  Pdm,t  Pvt+a  -pdmt)  Pvt+(l-Pvt)\ 

m- 1  '  ' 

Expanding  the  bracketed  expressions,  this  simplifies  to: 

M 

E  P>nmp-PsmJ  Pdml  pv] 

m=l  1 

This  quantity  is  represented  by  the  symbol  Pt  r 
Thus  the  prior  probability  of  the  event  is  written  as: 

/  \  Psi  Pdi  Pmi  Pv,  T 

f(eAis,z,o  =  n  -  '  _  '  '  n«,  (A5) 

i,eA|f,,0  l  pt  y=l  1 


Substituting  equations  (A3)  (A4)  and  (A5)  into  (A2)  gives: 

/  v  f,  (y,)  Ps,  Pdy  Pm,  Pv.  c  ( xr 

p(qa,u\z,o)  =  n  — _ —i  _ 1  ii  n  c  c,i 

g  (y  )  J  Pt  c=l  \  nc~nc,fi 

Cl  1 


YIpcc”c a-Pcf‘~n‘  UiM  Uptt  p(s\z,o) 

C=1  ;=1  / 


The  bracketed  expression  in  (A6)  is  constant  for  all  events  and  is  given  the  symbol  K.  The  noise 

density  term  in  the  initial  product  may  be  rewritten  as  Sc‘  — .  The  Pc  in  the  denominator  of 

Pc  ci 

ci 

this  expression  then  cancels  with  the  Pcc  c  term  in  the  second  product.  Similarly,  the  Ps}  term 

in  the  first  product  is  cancelled  by  the  denominator  of  the  target  density  function  fx  (y;) .  Thus 
equation  (A6)  simplifies  to: 

/  v  „  fx  (y,)  Pd,  Pm,  Pv.  c  (  a t  -n  \ 

p{qa'U\z,o\  -  k  n  —  1  1  1  n  c  cA  (A7) 

v  ’  VAI^o  gc{y)IPth  i\\nc~nc,  a)  (A7) 
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define  =  \ 


fm,&)  Pdm,tPnlm,tPVt 

Zcfy?  i 

i 


fe[l,T],  me[l,A4] 

f=0,  m=0 


then  we  may  write: 


/  v  JL  c  (N  -n  .)! 

p(eA,Lr|z,o  =  K  nw»,.ip  n  c  cA 


/=i 


c-1  <«c-”c,A)! 


(A8) 


Substituting  equation  (A8)  into  (Al)  then  gives  the  posterior  probability  of  the  joint  event 
defined  by  the  index  A. 

We  now  proceed  to  derive  an  expression  for  the  marginal  event  probabilities. 

Recall  that  for  ie[l,I] 


so. 


0.  =  U  0. 

i,m  A 

A  |A.k=  (l,m) 


Pr  (e,  „|  U.O,z)  =  Pr 


U  0A 

li,0,z 

k  A  |  XK=  (1  ,m) 

E  P>{eA\u,o,z) 


A  |AK=  (1  ,m) 


(A9) 


^  IL  c  (AT  -n  .)! 

e  nv'jy,-1,)  n 


A  |A.k=  (l,m) 


t  Wc-nc)!  (nc-ncA)\ 


and  similarly. 


Pr 


M 


U  0n  U0  1 

^  0  ,m  -1  ,n 

\  m= 1 


/r 


e  nt/y/^  n 


c  (M  -w  .)! 

c  c,  A 


A|^1V;  /=! 


^  Wc-«c>!  K-\a>! 


(A10) 


If  we  define  bi  m  such  that 
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pfej  u,o,z) 


M  I 


£  £  bj,n 
m=l 


then  a  solution  for  b,m  for  if  [1,1]  as: 


e  nw;  n-^- 

A|\=d,m)  /=1  '  C,1  ( Nc-nc)\  (nc-ncA)\ 


(All) 


The  summation  above  is  expanded  using  the  track  lists  described  in  section  3.1  to 


give: 


b. 


i,m 


E E -E  E  E  E  -EE  ft  W  ft ~ 

fi  mi  Vi  "Vi  mKt j  tIT  mIT  (=1  '  c=l  (N  -« 


^c-«c,a)! 


c  V>!  ^c-”c,a)! 


(A12) 


For  compactness  the  summation  ranges  in  ( A12)  have  been  suppressed.  The  summation  ranges 
implicitly  demonstrate  the  nested  summation  structure  since  the  range  for  each  measurement 
depends  on  the  assignments  of  all  the  preceding  measurements.  These  summation  ranges  are 
defined  as  follows: 


w 


££  -  £ 

'/  mi  frV0) 


W 

£ 

mrvi4X) 


Now,  it  is  possible  to  factorise  (A12)  since  the  depend  only  on  the  summations  prior 

to  /.  That  is,  m,  and  t,  are  independent  of  all  mp  tp  for  p>l.  Also,  due  to  the  marginal  event  of 
interest,  rn^and  tK  are  always  m  and  1  respectively  independent  of  all  other  m„  t„.  So  we  may 
rewrite  (A12)  as:  p  p  1 


>E E E  EEt. 


•"EEV.W  n  Wc'"‘a>! 

t,r  mIT 


(A13) 


C=1  (Nc-nc)l  (nc-ncA)l 


Now,  since  mp  tp  are  independent  of  all  model  associations  m,  for  other  tracks,  we  may  collect 
together  the  ij/m  as: 


Substituting  this  gives  the  final  expression  for  bjm  as: 
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Km  =  My1d1).  •  -E  ^(yK-i^K-i) 


K-l 

c 


•••E'Ky/rV  n 


Sc+l 

^c-«c,a)! 


(A14) 


C=1  Wc-”c)!  (”C-nc,A>! 


Similarly,  for  the  other  events  we  obtain: 

c  (N  -n  )! 

^=P''(0O,J0-UU-)E^1)  •••Ew  fl  c  C,A 

trr 


c- 1  GVWe>!  <V”C,A)! 


Pm,  PvAl-Pd,  Ps,  )  _  _  (JV  )[ 

_  i*  i  n  C  C'A 


(A15) 


Pt, 


c- 1  CWc-«P!  (Hc-na)! 


and 


__  c  (JV  -n  )! 

&-l,m  =  E  Wi}  -E  n  T^T  M 

t,  trr  c= i  (N„-hj!  J! 


Pm,  (1  -  Pr;,)  „  k  _c_ 

=  --  ,m  - E  t(yvfx)  -E  t(y/T,f/T)  n 

Pt  ~  fl  f/T  C-l 


'c  'V  v"c  c, A 


^c-\a>! 


(A16) 


(Nc-nc>!  <V"c,A>! 


As  previously,  the  summation  ranges  have  been  suppressed  for  compactness.  It  is  important 
to  note  that  these  equations  are  nested  summations  and  not  simply  the  product  of  sums.  That 
is: 

E  My  yd  E  'Ky^E  Myyt3)(. .  E  Myn,tir)  n  ---■-■)•  •  • 

*i  V  *2  ^3  trr  c=i  (Nc-nc)l  (nc-ncA)l' 


rather  than 


(E  'Ky^XE  Myrt2)fE  Myyt3)\  •  -(E  MyIT,tIT))H  (Nc  ”c,a)' 


=i  c Nc-«c)!  («c-nC/A)! 
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